sce <- readRDS(file.path(path_to_data, "sce_original_RPs_MPs.rds"))
sce_b <- subset_sce(sce, "activation", "baseline")
sce_a <- subset_sce(sce, "activation", "stimulated")
de_res_b <- readRDS(file.path(path_to_data, "de_results_baseline_original.rds"))
de_res_a <- readRDS(file.path(path_to_data, "de_results_stimulated_original.rds"))
features <- "state"
activation <- "baseline"
plot_marker_distributions(sce_b, features, activation)
## [[1]]
##
## [[2]]
##
## [[3]]
features <- "type"
activation <- "baseline"
plot_marker_distributions(sce_b, features, activation)
## [[1]]
##
## [[2]]
##
## [[3]]
features <- "state"
activation <- "stimulated"
plot_marker_distributions(sce_a, features, activation)
## [[1]]
##
## [[2]]
##
## [[3]]
features <- "type"
activation <- "stimulated"
plot_marker_distributions(sce_a, features, activation)
## [[1]]
##
## [[2]]
##
## [[3]]
perform_and_plot_DR(sce_b)
## [[1]]
##
## [[2]]
##
## [[3]]
##
## [[4]]
perform_and_plot_DR(sce_b)
## [[1]]